CONDITIONING OF QUADRATIC HARNESSES 



WLODEK BRYC AND JACEK WESOLOWSKI 



Abstract. We describe quadratic harnesses that arise through the double 
sided conditioning (bridges) of an already known quadratic harness and we 
characterize quadratic harnesses that arise by this construction from Levy 
processes. We also analyze a construction that produces quadratic harnesses 
by "gluing together" two conditionally-independent g-Meixner processes. Our 
main tool is a deterministic time and space transformation which we use to 
classify a class of harnesses with quadratic conditional moments. 



1. Introduction and main results 

1.1. Quadratic harness property. Throughout the paper T = {J-s,t) is a fam- 
ily of sigma fields with s < t from a nonempty open (generalized) interval T = 
(Tq, Ti) C (—00, 00) such that J-g t C Jv.u for r, s, t, u G T with r < s < t < u. We 
include Tq = —00 or Ti = cxd among the possible choices for the end-points of T. 

An integrable stochastic process X = {Xt : t G T} is called a harness [131 [ISl HO] 
on T with respect to T if Xu is J-'s^t measurable if u > < or if u < s, and for any 
s,t,u £ T with s < t < u, 

11 — / / — e 

(1.1) E(Xt|J-,,,0 - X, + -Xu. 

u — s u — s 

All integrable Levy processes are harnesses with respect to their natural filtration 
( |151 (2.8)]); additional examples appear in references on quadratic harnesses that 
are mentioned after Definition 11.11 

For a square-integrable process, a natural second-order extension of (jl.ll) is the 
requirement that Wa,T(Xt\J-s,u) is a quadratic function of Xs,Xu. It turns out 
that this assumption is much more restrictive than (|1.1|) alone. Wesolowski fl9j 
determined that there are only five such Levy processes. Under certain additional 
assumptions [31 Theorem 2.2] asserts that there exist numerical constants r],d eR 
(T, T > and 7 < 1 + 2^/ctt such that for all s < t < u, 

(1.2) Var[X,| J-,,„] = ( 1 + ^-^^ ~ ' ^ 



u ~ s u — s 

iuXs-sXu)^ {Xu-Xs)^ AXu-X,){uX,-sXu 
+r—, y2 (1-7)- 



{u — s)^ (u — s)^ (it — sy 
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where 



(1.3) 



Ft,. 



(u-t)(t-s) 



u(l + sa) + r — S7 



In this paper, we adopt these formulas as the definition, but fohowing W, Corollary 
5.1] for finite intervals T we allow negative values of parameters a, r, and we allow 
any 7 > ~1- 

Definition 1.1. Let T = {To,Ti) C (0,oo). We will say that a square-integrable 
stochastic process X = {Xt)ti£T is a quadratic harness on T with respect to (J>,s), 
if Xu is J^s,t measurable if u > i or if u < s, 



and equations (jl.ip . and (|1.2|) hold. When we want to indicate the parameters, 
we shall write X e Qi7(77, 0; ct, r; 7). (Formula (|1.3p then follows by taking the 
expected value of (|1.2I) . see Proposition 13.21 ) 

Examples of quadratic harnesses on (0, 00) are five Levy processes with quadratic 
conditional variances from [19j . Other examples include the classical versions of 
certain free Levy processes ([5j Theorem 4.3]), classical versions of (7-Brownian 
motion (fF, Theorem 4.1]), bi-Poisson process [31 [6l [7] and Markov processes with 
Askey- Wilson laws (II4 Theorem 1.1]). 

Remark 1.1 (Caution). While the use of general time interval (Tq, Ti) is convenient 
for the purpose of studying transformations, the reader should be aware that most 
of the fundamental results (uniqueness, integrability, orthogonal polynomial mar- 
tingales) valid for (0,(X)) do not extend automatically to quadratic harnesses on 
finite intervals, and may require additional assumptions on the interval and on the 
process. For example, from the of proof of [21 Theorem 4.1] one can deduce that the 
moments of a quadratic harness on (0, T) with —1 < 7 < 1 — l^far are determined 
uniquely if the conditional moments of Xt and X"^ with respect to past cr-field Ts 
for s <t are linear and quadratic functions of Xs respectively. This result fails for 
intervals bounded away from 0, see Example 11.21 or when the assumption on one 
sided conditioning is dropped, see Example ll.il 

1.2. Results. Our main result shows how to transform a generic harness with 
quadratic conditional variances and with a product covariance into a quadratic 
harness. 

Theorem 1.1. Lei H. be a harness (|l.ll) with respect to the family T on an interval 
{Tq,Ti) C II with mean 



(1.4) 



E{Xt) = 0, nXtX,) 



min{<, s}. 



E{Xt) 



a + pt 



and with covariance 



(1.5) Cov{Xs,Xt) ^ {as + b){ct + d), s<t, 

such that ad — bc> and {at + b){ct + d) > on [Tq, Ti). Suppose that 




) 
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where Ft^s.u is non-random, and 

(1.7) x-=X + ari + ei3 + (ja^ +tI3^ +pal3>{). 

Let il}{t) = {dt — h)/{a — ct). Then stochastic process 



(1.8) 



a — ct 
ad — he 



X. 



is a quadratic harness in QH(r]' ,9']a' ,t']^') on the interval ^ cTo+d ' ? 
(0,oo), and has parameters 



Ti+b 
cTi+d 



C 



(1.9) 
(1.10) 
(1.11) 
(1.12) 

(1.13) 



a 

T' 
1' 



{d{ri + /3p + 2aa) + c(i 
{b{r] + Pp + 2aa) + a{: 
{tc^ + dpc + d^cr)/x 
[ra^ + hpa + 6V)/x 
bcp + adp + 2hda - 



1 



ap - 
ap 



2acT 



X 



Theorem 1.1 was motivated by a re-parametrizations of orthogonality measures 
of the related families of orthogonal polynomials when ct = 0, which was discovered 
by R. Szwarc; this observation is presented in an unpublished manuscript [5]. A 
version of Theorem 11.11 is implicit in the construction of quadratic harnesses from 
Markov processes based on the Askey- Wilson integral, see 1, Formula (2.28) and 
Section 6.2]. 

For reference, we state an elementary special case: if X G QH{ri,9;a,T;^) and 
Zt — aXi/^2 with a ^ 0, then 

(1.14) Z e QH{ri/a, aO; aja^.a^T; 7). 

A related transformation Yt — tX^jt produces Y e QH{9,ri;T, a;'^), i.e. entries 
within pairs {r],0) and (ct, r) are swapped. In fact, these two elementary facts can 
be used to re- write formulas (|1.9fll.l3|) into other "equivalent forms" . 

We also remark that the transformation used in Theorem 1 1.1 1 is reversible, so X 
can be "represented" in terms of the quadratic harness Y e QH{r]', 9'; ct', r'; 7') as 

Xt = {ct + d)Y(^at+b)/{ct+d) + a + (3t. 

Next, we apply Theorem 11.11 to analyze/ classify processes that arise through 
double-sided conditioning. For a (random or deterministic) real function X of real 
parameter t, denote 

rA.,n(x)i 



A,,JX 



with 
(1.15) 



A,,4X) 



X,, — Xx 



A,,„(X) 
and As,u{X) : 



uX, - sX„ 



In the sequel, if X is clear from the context, to shorten the notation we will 



sometimes write ^ instead of ^ (X) 



The right hand side of ()1.2|) is a quadratic polynomial in two real variables which 
we will write as 



(1.16) i^(A, J = 1 + ?/A,,,„ + ^^A,,, + aAi 



rA 



(1 - 7)As,„As,„ 
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It will be convenient to write if (a, b) instead of K 



Theorem 1.2. Suppose {Xt)t£{To.Ti) G QH{7],9;a,T;j). Choose R < V in {Tq,Ti). 
Then V{1 + Ra) + t - Rj > and 



M = M{Xr,Xv) - ^K{^jiy{iq)/^V{Ra + 1) + r - i?7 
is non-zero with probability one. Define 



(1-17) Yt = jy^^^^ ((1 + - tXviV - Xr) 



for t > 0. Then for < s < t < oo, conditional version of (jl.4p holds: 

E{Yt\TRy) = 0;EiYsYt\TR,y) = s. 
With Gs^u = T(s+R)/(i+s/v).(u+Ry(i+u/v), condition (fLTj) holds: for s < t < u, 

EiYtlGs.u) = ^Y, + *—^Y^ 
u — s u — s 



and condition (|1.2p holds: for s < t < u, 
(1.18) Var[r,|^,,„] = ( 1 + " 



u — s u — s 

2 ^v- _ ^2 



(u — s)2 (u — s)2 (u — s)^ 

with J-'r^ -measurable parameters: 

(1-19) 

^_ -0 + Vv- 2TAR,y{X) + 2aVARy{X) - (1 - 7)(V^Afl,v(X) - Afl^v(X)) 
^ VVM{V{Ra + 1) + T - i?7) 

(1-20) 

^ Vv{e-Rr^ + 2TARyiX) - 2RaARyOL) - (1 - j){ARyiX) - i?A,j,v(X))) 
^ ~ M(y{Ra + 1) + T - i?7) 

and wzi/i non-random 

aV^ + {l-j)V + T 



(1.21) 



(1.22) 
(1.23) 



y(F(i?(T + 1) + T - i?7) 

cri?2 + (1 - 7)i? + r 



y(i?o- + l)+T-i?7' 
^7 - i?(y(T + 1) - T (F-i?)(7 + l) 



V{Rcr + 1) + T - i?7 i?cry + y - i?7 + r 

(Constant Ft^s^u is then given by formula (|1.3p with a,T,j replacing a,T,^.) 

We now specialize this same result to Markov processes. Notice that quadratic 
harnesses have cadlag versions, so regular versions of conditional distributions ex- 
ist. If Z is a quadratic harness with respect to {J^s,t} on (To,Ti) and R < V 
are in (TqjTi), then conditionally on Z{V),Z{R), the process still has quadratic 
conditional variances with respect to J^s,u when s,m e {V,R), and of course the 
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conditioned process takes deterministic values at the endpoints V and R. Further- 
more, the conditional variance ()1.2p is positive with probability one. Thus there is 
a set of probability one of pairs zjf e supp(Z(i?)), zy € supp{Z {V)) such that the 
laws 

(1.24) TTtiU) := Pr(Z(<) G U\Z{R) = zr, Z{V) = zy) 

are well defined for all t G (V", i?) and that there are Borel sets Us of TTs-measure 
one such that 

(1.25) Ps.t{x, U) Pr(Z(t) e U\Z{s) = x, Z{V) = zy) 

are well defined for all s < i in (V, i?) and all x G Us- It is known that (|1.24p and 
(|1.25|) determine a Markov process if Z is Markov, or more generally, if Z is the so 
called reciprocal process, see [171 Theorem 4.1]. 

The following is the corresponding re-statement of Theorem 11.21 

Theorem 1.2'. Let Z = {Zt)te{To,Ti) e QH{q,9;a,T;^) he Markov. Fix R < V 
in {T{),Ti) and z^ G supp(Zfl), zy G supp(Zv) such that with A^y :— [zy — 
zr)/{V - R) and Any := (Fzi^ - Rzy)/{V - R), we have K{ARy, A^y) > and 
such that (fTMl and (fOSj) are well defined. Let X ^ {X{t) : R < t < V} be 
the conditional process, i.e. the Markov process with univariate laws (|1.24p and 
transition probabilities (|1.25p . 

Then V{1 + Ra) + t - Rj > and 

= K{ARy, ARy)/{V{l + Rcr)+T- Rj) > 0. 
Process X is related by a deterministic transformation 

(1-26) Y, = (^T^ly^ {{t/V + l)X(,+«)/(,/y+i) ^ (t^ + z^)) 

(compare p.l7p ) to i/ie quadratic harness Y G QH{rj,9;(7,T;^) on {0,oo) with 
parameters given by U.2(MT7S3\) with A^y(X) = [AflVjAi^y]-^ which is now a 
non-random vector. 

Remark 1.2. Analogous results with essentially the same proof hold for one sided 
conditioning on (0, V) and on (i?,oo). Formulas for parameters in these two cases 
correspond to Theorem 11.21 after taking the limit as i? — > or as F — > oo, 
respectively, with conventions that limy_>ooAfly — 0, Huir^q Ary ~ zyjV, 
limy^oo ARy = zr, limfl^o A^jy = 0. More specifically, suppose Z is a quadratic 
harness on (0, oo) in QH{ri, 9; cr, r; 7). 

(i) Let X = {X{t) : < t < V} be the conditional process arising by condi- 
tioning Z with respect to Zy — zy. Assume that 

Then 

(1-27) Y, = (^(1 + i/V)X,/(i+,/y) - Izy) 
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is a quadratic harness with parameters 
6i + 2rzy/l/ 



(1.28) 0Y = 

(1.29) 7JY = 

(1.30) TY = 

(1.31) O-y = 

(1.32) 71- = 



-e + Vrj- 2tzv/V - (1 - 7)zy 
Ft 



V + t' 

V{V + t) 

— T 



(ii) Let X — {X{t) : R < t < oo} be the conditional process arising by 
conditioning Z with respect to = zr- Assume that 

= {l + Ra){l + rjZR + az\) > 0. 

For < i < cx), let 

Yt = {l + Ra){Xt+R- zr)/k. 
Then Y is a quadratic harness with parameters 



(1.33) Oy = 

(1.34) ryy = 

(1.35) Ty = 

(1.36) CTy = 

(1.37) 7y = 



0-Rrj- 2RazR - (1 - j)zr 

K 

<tR^ + (1 -7)i? + T 
l + Ra ' 



l + Ra 
7 — Ra 



l + Ra 

Conditioning ahows us to identify new quadratic harnesses, as seen from the 
foUowing. 

Theorem 1.3. Suppose that cr, r > are such that ut < 1, ^ — 1 — 2y/aT and 77, 9 
are real numbers such that ^frr] + ^/ct^ = 0. Then there exists a square-integrable 
Markov process (^t)tG(o,oo) such that (jl.ll) . (jl.2p . and (|1.4p hold. 

The range of parameters in Theorem I 1 . BI corresponds formahy to g = 1 in [T|, and 
perhaps could be obtained from the processes in that paper by a limiting procedure. 

Proofs of Theorems 11.11 and ll.2'l are in Section 13.21 Theorem 11.31 is proved in 
Section 1121 



1.3. Examples. Here we use Theorem 1 1.1 1 to give examples that illustrate lack of 
uniqueness for quadratic harnesses on finite intervals. 

Example 1.1. Let {Wt) be the Wiener process and ^ be a centered random variable 
independent of W with E^^ = v"^ . Let 

Xt^Wt+ ^t, t > 0. 
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Then E{Xt) = and Cov{Xs, Xt) — s(l + v'^t). Furthermore, Xt is a harness with 
respect to its natural sigma fields, and 

WsiT{Xt\Ts,t) = YMWtlWs, Wu) = Ft^s,u. 
So from Theorem 11.11 (or by direct calculation) we see that 

is a quadratic harness on (0, 1/w^) with respect to natural cr-fields and has param- 
eters r] = 9 = a = T = Q,j — 1. 

Example 1.2. This example is a time-inversion of Example 11.11 Let (Wt) be 
the Wiener process and ^ be a centered random variable independent of W with 
E^^ = v^. Then, with respect to natural a- fields, Xt = Wt-v2 + ^ is a quadratic 
harness on (u^, oo) with parameters r] = = a = T = 0,j = l. It is also easy to 
see that Xt and Xf — t are martingales with respect to the past cr-field. 

Next, we give a simple example of a quadratic harnesses with 7 > 1 and err > 1; 
such examples are interesting because most of the general theory developed in [3] 
does not apply. 

Example 1.3. Suppose {Gt)t>o is a gamma process with parameters (1, 1) (see 
Example 14.11 below) . Let ^ be an independent random variable with mean E{^) — 
^ > and E^'^ = Let 

Xt = ^Gt 

Then E{Xt) = (3t. From 

CoviX,,Xt) =EiCoviX,,Xt\0) +CoviEiX.,\OMXt\0) 
we see that for s < t, Cov{Xs, Xt) ~ s{v'^t + j3). Let J^s,u be the natural cr-fields 
associated with Xt. Consider the auxiliary cr-fields J^s,u generated by ^ and {Gt : 
t e (0,s] U [u,^)}. Then E(Xt|J^,,„) = ^ (^G, + ^G„) = + so 

E{Xt\Ts,u) = ^X, + ^Xu. Similarly, using ^ with a = 1 we get 

Wai{Xt\Ts,u) = e^s.TiGt\G,,Gu) = {n-t){t -s) _ ^^^^ 



S + 1)(U - S)2 

{u-t)it-s) 
{u — s + 1)(m — s)2 



{Xu — Xs 



so 



Var(X,|J-,^„) = ^-A!i-^1^^^-^(X„ - X,)\ 



(w — s + l)(w — s) 
From Theorem 11.11 applied with 

a ^ V, 6 = 0, c = V, d = f3/v 

we see that 

Zt = v{l - t)X Jit - —t 

i.^(l-t) V 

is a quadratic harness on (0, 1) with parameters 

T] = e = 2v/l3, a = T = 7=1 + 2y/^. 

In particular, 7=1 + and ctt = v^/P^ = (^'(C^)) > 1 can be 

arbitrarily large. 



8 



WLODEK BRYC AND JACEK WESOLOWSKI 



Of course, the distribution of is arbitrary so the moments of Zt are not deter- 
mined uniquely and may fail to exist. 



An example of quadratic harness on (0, oo) with 7 = 1 + is a bi-Pascal 

process [16j . which is a quadratic harness on (0,oo) with arbitrary cr = r > and 
arbitrary rj = 9 > l^fr. A related quadratic harness on a finite interval (0, T) arises 
as a transformation (|1.8p of the generalized Waring process introduced in [S] and 
[S], see also [H] and [H]. 

2. Matrix notation 
For calculations, it will be convenient to parametrize time as a subset of a pro- 

We rewrite in vector form as 



jective plane, i.e. using t 



(2.1) nM^s.u) - (t, A,^JX)> , 

where (a, b) — a'^b, and the components of A ^ „ (X) are defined by (|1.15[l . 

It follows from (jl.ip that admissible expectations of a harness X are afhne in 
i.e., 

(2.2) E(X0 = (t,M), ter, 

where 



(2.3) 



Moreover, if X is a square integrable harness then by [3, Proposition 2.1] the ad- 
missible covariances are of the form 



(2.4) 
where 



Cov(X„ Xt) = (s, Et) , s,t(^T, s<t 



Co 

C2 



Cl 

C3 



and s — 



Throughout this paper, letters s,i, u G T are reserved to denote time, and s, t, 
and also u = [u, 1]^ have this special meaning also when used with subscripts or 
primed. We also use the convention that s < t < u. 

Note that under our convention s < t so T, is not a symmetric matrix; for 

example, covariance s At is represented by matrix S = 



Formula (|1.2|) in matrix form can be written as 

(2.5) Var(Xt|J-,,„) = (l + {0,A^J -t 
where 

(2.6) t 







(A FA 



)) . 





, r = 


r -1" 










7 a 



Here 77, 0, cr, r, 7 are constants independent of s, t, u. 

Of course matrix F is determined only up to Fi_2 + r2,i- The usual choice of 
symmetric F is in fact inconvenient, see Proposition 13.21 The choice made in (12. 6p 
matches the notation we used in previous papers: after substituting q for 7, the 
resulting parametrization of the conditional variance is identical to O (2.14)]. 

The non-random constant Ft^s,u is determined uniquely by taking the average 
of both sides of (|2.5p . According to [3l (2.15)], with the choice of F as in p.6p . 
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formula ()1.3p holds. We re- write formula (|1.3p in matrix notation using a special 
matrix 



(2.7) 



J 



1 

-1 



It is easy to see that = — /, = — J. For future reference we state also two 
less obvious properties: for A E GL2(B,), 

(2.8) A^JA = det(A) J and j'^AJ = det{A){A-'^f . 



Formula (II. 3p can now be written as 



(2.9) 



(t , Ju) (s, Jt) 

(s, j^r ju) 



This formula makes sense for any 2x2 matrix F as long as the denominator is 
non-zero. 

Finally, we note for future reference that the conditional covariance for quadratic 
harnesses also takes a simple form: with s < ti < t2 < u, 



(2.10) 



where K{a) = 1 + {d_, a) + (a, Fa) is the quadratic polynomial from (|1.16|) and (|2.5p . 
A quick way to see this is to notice that implies 



U ~ t-2 
U — tl 



Var(XtJJ-,,, 



(ti, Ju) 



Var(XtJJ-,,„). 



3. Deterministic time and space transformation 
With a non-degenerate afhnc function / : — > R^, written in matrix notation 

as 

(3.1) f{x,y) = [x,y]A+[mi,m2] , 

we associate Mobius transform ip{t) — {at + h)/{ct + d) generated by 



A = A 



f 



a b 
c d 



e GL2{R). 



If X is a stochastic process on an open interval T C R, and T lies in the range 
of (fi, we define a deterministic transformation X-' of the stochastic process X as 
the process Y = X^ on 5 = (p~^{T) such that 



(3.2) 



Yt = {ct -f d)X, 



(t,m), t e S. 



We note that (p is increasing on S if det(A) > and it is decreasing otherwise. 
Our interest in this transformation comes from the fact that special cases ap- 
peared as ad hoc tricks in constructions of quadratic harnesses in [BJ [1] [16] . 

A calculation verifies that as long as the time domains of the processes match, 



g°f 



(X^)5 = X 

This allows us to build more complicated transformations in simple steps, and gives 
us flexibility to consider either Y = X-' or X = Y-^ as needed. 
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For non-degenerate afRne / with Mobius transform : {To,Ti) 
s < t in T , the transformed cr-fields are 



[T^Ti) and 



(3.3) 



j-f _ I ^^(s),ip(u) 



if det(^) > 0, 
if det(^) < 0. 



It is clear that if X has hnear regressions and quadratic conditional variances with 
respect to J-s,u then has linear regressions and quadratic conditional variances 
with respect to The following technical result describes how the parameters 

transform in the setting slightly more general than Theorem 11.11 

Proposition 3.1. Let 'X. be a harness (jl.ip with respect to the family T on an 
interval {Tq,Ti) with the first two moments given by (j2.2p and (j2.4p . Suppose 

(3.4) Var(X4|J-,,„) = {x+{i,As,u) + (A,,„,rA,^J) , 

with non-random Ft^s,u, X € K-; £ G ^ j '"^'^ arbitrary 2x2 matrix T . Let f be 

a non-degenerate affine function with matrix A — ^ ^ , corresponding Mobius 

transform ip and shift m. If ip is well defined on the entire interval {To,Ti) then 
X :— X-^ given by 

Xt = {ct + d)X^at-i-b)/{ct+d) + m-it + 7712 
also satisfies (P^ . ([23) with 



(3.5) 

(3.6) S 
and pTT]) and (I33) with 

(3.7) f = 

(3.8) 
and 
(3.9) 



A^Y.^A 



m. 



A-^Y{A-^Y 
[A-^fVA-^ 



if det(A) > 0, 
%f det(A) < . 

if det(A) > 0, 
if det(A) < 0. 



A-^d_-ij: ^V^)vci. 



X = X - (^,m) - (m,rm), 

We remark that in the most interesting case of "product covariance" 
(3.10) Cov(X^,Xt) = (es + (?)(# + i/'), s<t, 

transformation p.6p preserves this product form. In fact, in this case 

but if we collect the coefficients of the covariance (|3.10p into another matrix 

e = P ^1 

then a calculation based on p.6p shows that for det ^ > the covariance of X 
corresponds to 9 = 6^. 
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We postpone the proof of Proposition 13.11 until Section 13.11 so that we can first 
clarify the role of non-random constant Ft^s,u- The main point is that in the non- 
degenerate case with ci > C2, this constant is determined uniquely by taking the 
average of both sides of (|3.4I) . Furthermore, Ft^s,u is often given by formula (|2.9I) . 
i.e., (|L3l) . 



Proposition 3.2. Suppose a harness X has mean ()2.2|) and non-degenerate co- 
variance (|2.4p with ci > C2 . //X has quadratic conditional variance (j3.4p and the 
off- diagonal entries of matrix T are chosen so that 

(3.11) x + (^,A^> + (Ai,r/i)+tr(rS^) = 0, 

then (s, J-^rju) = u(l + sa) + r — 57 7^ and Ft.s,u is given by formula (|2.9I) . 
Moreover, transformation formulas in Proposition \3.1\ preserve 



Formulas (jl.2p and (|2.6p illustrate the choice of such F. 
3.1. Proof of Propositions [311 and [3721 

Lemma 3.3. Let f be a non- degenerate affine function (13. ip with Mobius transform 
Lp. If s' — ip{s), u' — (p(u), then 

(3.12) A,,JX-^)=A^A,,,„,(X)+m. 

Proof. Let g{x,y) = denote the linear part of /. Since Aj, ^,(a) = m on a 

linear function a{t) := (t,m), and X-^ (t) = X^ (i) -I- (t,m), we have „(X-^) = 
Ag ^(X^) -I- m. Since Xs' = X^(s)/(cs -I- d), and from the matrix form of (|1.15p we 
have 

(3.13) A.,(X).^(^*^-^^^) 



(t , Js) 



we get 



A ex] — T cu+g — cs+a — _ j 

— "''"'^ ' ^ (u',Js') " ((cu + d)u', J(cs-Fd)s') 

Noting that 

(3.14) (cs + d)s' = As, 

and using (j2.8p we get 



(u, A^J^s) ^ ' (u,A^JAs) 

Thus A,,„(X/) = A,,„(Xff) + m = A^A,,,„, (X) -f m. □ 

Proof of Provosition lS. 11 Throughout the proof we write t' = ip(t). If if is increas- 
ing, by (j2.ip and the definition of X-'' we have 

(3.15) E(X-^(i)|J-/^J = {ct + d)E{Xt>\Ts',u') + (t,m) 

= (rf + d)(t',A,,,„,(X)) + (t,m). 
By ((3111) and ([3l2l) we get 

E{Xf{t)\:FlJ = {At, (A-^ (A,_„(X/) - m)) + (t,m) = (t, A,,JX/)) . 
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Thus the condition holds true and is a harness. Similarly, one can verify 
that (II. ip holds when (p is a decreasing function. 
We use (|3.14[) to compute the mean of X-'' 

E(X-^(<)) = {ct + d){t',ti) + (t,m) = {At,fi) + (t,m) = (t,m + A^/i) , 



and (|3.5I) follows. 

To find the covariance we again use p.l4p and the fact that Cov {Xs' , ) is 
either (s', Et') or (t', Ss') = (s', E^t') depending whether s' < t' (case det(^) > 0) 
or s' > t' (case det(A) < 0). For example, if det(A) > then 

Cov(X-'^(s),X^(t)) = {cs + d){ct + d)CoviXs,,Xt') 

= {cs + d){ct + d)(s', Et') = {As, SAt) = (s, A^EAt) , 

and thus p.6|) follows. 

Since is monotone on (TqjTi), 



(3.16) Var(Xt|J-/.J = 



■ (ct + dfYniiXf I J"/, „,) if det(A) > , 



{ct + dfWaiiXt'lTl, ,^,) if det(A) < . 



Consider the case det(A) < so that u' < t' < s' . Since A„ j, = Ah,a^ by p.l6p 
and Lemma [3.31 the conditional variance is 



(3.17) Yar{Xt\FU = {ct + df Ft, (x + {A'^gx, ^sA^^) ^ ^) 

+ (A,,JX/) - m, A-ir^(A-i)^(A,,JXO - m)) 

So p.l7p rewrites as 

Var(Xt|J-/j - (x - (A-i^x,m) + (A,, JX), A-^^^) 

+ (A,^„(X) - m, A-ir(A-i)^ (a,^JX) - m)) 
= -Ft',n',.'(x- (^"'^x,m) + (m,r^m) 
;A,,JX),A-i^- (f + f^)m\ + (A,,JX),f A,.JX) 



(s,Jt)(t,Ju) 1 , /^_(x),^-i,_(r + F)^ 



(s,j^rju) G 



+ (A,_JX),r^A,^JX)) 

Since the last term is invariant under transposition, we get p.8p and (|3.7p . The 
case det(A) > is handled similarly and the proof is omitted. □ 



The proof of Propistion 13.21 is based on the formula for the covariance matrix of 
vector Ag.u- 

Lemma 3.4. 

(3.18) Cov A, „ = ~ Jus^J^ + E^. 

u — s 
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Proof. From p.l3|) we get 

Cov A,,„ = E (a,_„aJ;„) - E (A,^ J E (a^„) 

[u — s)2 

Note that smee S"^Sw — u^Yj^s and su^ — us^ — {u — s)j'^, the numerator can 
be written as 

J{su^ - us^)i:us^j'^ - J{s%F - us^)lfsu^J^ = (w - s'^ - ifsu^).]'^ 

(recall that = — J and JJ^ ~ I). Further we write the above expression as 
(w-s)[(I]-S^)ws^+S^(us^-su^)]j'^ = {u-s){ci-C2)Jus^J^ + iu~s)'^T.^JJ^ 
and thus ([XTH]) follows. □ 
Proof of Proposition We first remark that formulas p.7p and p.6p imply that 

tr(ri]^) = tr(fE). 

Next, we note that 

X + H^Hl) + it, Tm) = X + il.D + % f^)- 

(This follows from a longer calculation based on the formulas from Proposition l3.9l ) 
Therefore, transformation formulas preserve p. lip . 

Next, we show that p.lip implies (|2.9p . This will be accomplished by computing 
the averages of both sides of p.4p . 

For any harness with covariance 



(3.19) E (Var(Xt| = IL^IIILJ) - c^). 

u — s 

To prove ((3T9)) . we use (|3J8)) . From ((2l1) we get 

EVar(Xt| J-,^„) = VarXt - Var(E(Xi | J-,^„)) = VarX* - fCoY{A^Jt 

On the other hand, with ^^(A^^^) = x+ {i,As,u) + (As.„,rA^ „), we have 

(3.20) E (if (A,_J) = tr(rS]^) + X(a.) + {c, - c,)^ll^^^^. 
To prove (|3.20p we note that E ^ — ^, so 

(3.21) EK{A,J^x + KO'^As.u+^Alu^As,u 

= X + + M^F/. + tr (FCov A,^J = if (/i) + tr (FCov A,_J . 
From (|3T8)) we get 

(3.22) tr (rCov A, ,J = ^^ tr (F Jus^J^) + tr (FS^) . 
Since 

tr (F Jus^J^) = tr {s^j'^TJu) = (s, J^FJu), 
([5:^ follows from ([5:^ and ([X^ . 
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Since E (Var{Xt\J-s,u)) — Ft,s,uIE (K(A^ ,^^)) , in the non-degenerate case ci > C2, 
and in particular in the non-degenerate product (jS.lOp with det > 0, formula 
p.l9p implies that E (K(A.^ „)) 7^ so Ft.s,u is uniquely determined. 

Moreover, while the quadratic part of K does not determine the entries of T 
uniquely, it is natural to choose a unique, perhaps non-symmetric, matrix F such 
that tr(rE^) + A'(^) = 0, i.e. ([XTT]) holds. Then Ft^s,u is given by formula ([^ 
which is just a matrix form of (jl.3p . 

□ 

Remark 3.1. It is interesting to remark that in the non-degenerate product case 
([XTUl) with dot e > 0, formula ([XTO)) takes the form 

. , XX (t-s)(u-t) 
E Var J-^ « = det 9 > 0. 

u — s 

Remark 3.2. If C3 > 0, ci > C2, and cqc^ > then the right hand side of (12. 4p 
indeed defines a positive definite function on T = (0, 00). To verify this, let a(s, t) = 
(s, St) for s < i and c(s, t) = a{sAt, sVt). For sq = < si < • • • < s„, we compute 

det[c(sj, Sj)]o<i.j<„ 

= (Cl - C2)"(s„ - Sn-l) . . . (S2 - Sl)si (c3(ci - C2) + (cqCs - cl)Sn) ■ 

Then [c(si, Sj)]i<i.j<„ is positive definite as a sub-matrix of the above. 

From the proof of Proposition 13.21 we can also read out that processes with 
Cl = C2 are degenerate in the sense that Xt is a linear combination of Xg, X^. 

3.2. Proofs of Theorems [iTl] and |1.21 



Proof of Theorem\r^ Let A = ^ 
We apply Proposition 13. II with f{[x,y]) ~ {[x,y 



d -b 
-c d 



so that its inverse is _B = 
- [a,/3])B to 



ac ad 
be bd 



, F 



T p 
P o- 



From the transformation formulas we get /i = 0, S = 



1 




and X 



X + ar] + 



op 



and 



T/32 + 2pal3 > by ([TTD. We also get 

b{r] + I3p + 2aa) + a{9 + ap + 2/3t) 
(i(77 + /3p + 2aa) + c{9 + ap + 2/3t) 



F = 



^2 ' bpa + b'^a \{bcp + adp + 2bda 

; [bep + adp + 26(i(T + 2acr) rc^ -I- rfpc + d^a 

The quadratic polynomial K remains unchanged if we replace F by 

p, _ To^ + bpa + 5^(7 —X 

X + + adp + 26(i(T 4- 2acT rc^ -|- dpc -|- d^a 

Rewriting the quadratic polynomial as 

K{^) - X + (I, x> + (x, r'x) - X fl + x) + (x, 4f'x) 



TO 



2aeT) 



X 



X 



we get the parameters as claimed. 



□ 



CONDITIONING OF QUADRATIC HARNESSES 



15 



Proof of Theorem ITFl To prove that V{1 + Rn) + t - Rj > we note that V i~> 
V{1 + Ra) + T — i?7 is a contmuous function on R which by Proposition 13 . 21 cannot 
cross zero on {R,oo). 

The remainder of the proof consists of computing the parameters of X so that 
we can apply Theorem ll.il First, from (|2.ip we see that the mean of X is given by 

^^^^ Next, from (piU)) we see that CoY{X,,Xt) = KP{V - t){s - R) with 



A 



RV 



3.23) M = , " > . 

^/V{1 + Ra)+T^R-f 

Finally, we observe that the conditional variance of X is given by the same formula 
as the conditional variance for Z. Therefore we can apply Theorem 11.11 with 

a = Arv, fi = Arv, X = 1, 

a = A/\/y, b = -RMW, c = -M/W, d = mW. 

(Other choices are also possible and lead to "equivalent" quadratic harnesses as 
in ()1.14p .) Under the above choice of a,b,c,d, formula (jl.Sp gives (|1.17|) . Since 
X — K{Arv , Arv) > 0, assumption ()1.7p holds, and the parameters of the resulting 
quadratic harness are as claimed. 

□ 

4. Conditioning of quadratic harnesses 

In this section we show how to apply previous results to analyze which qua- 
dratic harnesses can be constructed from already know quadratic harnesses either 
by conditioning. Our basic building blocks will be the q-Meixner processes: these 
are quadratic harnesses with rj = a = and 7 = (7 G [—1, 1], see [5J. In particular, 
the five classical Levy processes mentioned in the introduction, sometimes called 
Meixner processes, correspond to 7 = 1. 

Our main interest in such constructions stems from the fact that the construc- 
tions from the Askey- Wilson laws [I] yield only 7 G (— 1, 1 — 2y/aT); in particular 
"classical" quadratic harnesses with the boundary value of 7 = 1 — 2y^aT that ap- 
pear in |3j Proposition 4.4] need additional work and have been constructed only 
for particular choices of the parameters. Example 14.11 and Example 14.21 present two 
new quadratic harnesses with 7 = 1 — 2^faT^ and in Remark 14.11 we mention two 
more cases. 

4.1. Conditioning of Meixner processes. We can use Theorem ll.2l to recognize 
which quadratic harnesses could arise by conditioning from quadratic harnesses that 
are Levy processes. 

Proposition 4.1. Fix t > and 6* e R. Suppose Z e QH{0,9;0,t; 1) i.e., Z is a 
Meixner process. If ^ is a conditional bridge of Z then 



is in QH{riY, Oy; ay , ty'tIy) with parameters 

V-R-T T 

= V-R + r ^ = = V-R + r ^ 



16 wlodek bryc and jacek wesolowski 

6 + 2t/\rv 



(4.3) By = -VY 



^V-R + T^l + 0Arv + rA-j^y 



Proof. From Theorem 1 1 . 2' I we read out the transformation that leads to parameters 

, V-R-Tz , tV , T 
7 = T7 ' ^ 



V-R + t' V-R + t' V{V-R + t) 

e + 2TAR,v 



VV-R + T^1 + OArv + TAly 



, e + 2tArv 



VVW -R + + QArv + rA^^^ ■ 
Further transformation (|1.14l) leads to (|4.1I) and to parameters as claimed. □ 

Example 4.1 (Dirichlet process). For any uq^tq > with (JqTo < 1, there exists a 
quadratic harness Y (namely, a Dirichlet process) which has parameters ay = cto, 
TY = To, By = 2^/7^, r^Y = -2^/5^, and 7y = 1 - 2y/aQTQ. 

Indeed, consider the case of conditional harnesses obtained from a gamma process 
{Gt)t>o- Gamma process (Gt) is a non-negative two-parameter Levy process with 
parameters a, /3 > 0. The density of Gt given by 

r*x"*-ie-'^7r(a)l(o,oo)(x). 

As a Levy process, (Gt) is a harness with mean E(Gt) = ta/(3 and variance 
Var(Gt) = ta/P"^. It is also known (see (|4.6p below) that 



(4.4) Var(Gt|J-,,0 = ^— ^—^^-^^^^ (G„ - G,)^ 



{t~ s){u~t) 
[u - s)2((m - s)/a + 1) 

Then 

■Z^t = —Gat ~ Oit 

a 

is a quadratic harness in Qi/ (0, 2/a; 0, l/a^; 1) which by further transformation 
(|1.14p can be transformed into an element of QH{0, 2; 0, 1; 1). 

Instead of considering a conditional process of (Gt) we therefore consider a 
conditional process of Z G QH{0,2;0,1;1) and apply Proposition 14.11 Choose 
V — R = ^(ToTo — o-qTo so that 1/{1 + V — R) = ^o-qTo- From (14. 2p we have 



Transformation (I1.14[) with = ^ ctq /tq gives 

Ty = ^JaoTo/a^ — tq, CTy = a^V'^'o'^o = ctq. 
By ((431) with any A := Arv > 

2(1 + A) 2 , 

W1 + 2A + A2 - i? + 1 

so 0y = 2y^. Similarly, rjY/a — —2^aoTo, so 77^ = —2^/ao. 

Since conditional processes of the gamma process are Dirichlet processes, the 
same conclusion can be obtained directly by a fairly natural reparametrization 
()4.7p without invoking explicitly any of the transformations. Let ai, . . . , a„, a„+i be 
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positive numbers. A Dirichlet distribution D„(ai, . . . , a„, a„+i) is defined through 
its density 

f{xi,...,Xn) = ""^^ n ^^"'"^ h l;7„(a;i,...,a;„) , 

where ?7„ = {(xi, . . . , a;„) e (0, cx))" : X]r=i < ^i- ^ stochastic process X = 
(Xt)tg[o,y] is called a Dirichlet process if there exists a finite nonzero measure /i on 
[0, V\ such that for any n and any < ii < . . . < i„ < F the distribution of the 
vector of increments {Xt^ , Xt^ — Xt^ , . . . , Xt^ — Xt^_-^) is Dirichlet 

Dn{fi{[0,h]),fi{[ti,t2]),...,n{[tn-l,tn]),f^{[tn,V])). 

This is one of the basic objects of non-parametric Bayesian statistics - see [TUJ [T^ . 
Let fi — cX, where A is a Lebcsguc measure on [0, V] and c=l/a>Oisa number. 
Recall that the beta distribution, Bi{a, b), is defined by the density 



and if X ^ Bj{a, b) then 
(4.5) E{X) = — and Var(X) 



a + 6' ^ ' (a + 6)2(a + 6+l) ' 

Since Xf has the beta distribution Bj{ct,c{V — t)) the formulas (|4.5p give 

E(Xt) = - and Cov(X„Xt)- ""^^ ^ 



Note that to compute K{XsXt) it is convenient to use the classical fact, that Xg/Xt 
and Xt are independent and X^/Xt is a beta Bj(cs,c{t — s)) random variable. 
Note also that X is Markov process with transition distribution defined by the fact 
that {Xt — Xs)/{1 ~ Xs) and Xg are independent, and {Xt ~ Xg)/ {1 ~ Xg) is beta 
Bj{c{t — s), c{V — t)). It is also known that 

Bl{c{t-s),c{u-t)) 

and {Xt — Xg)/{Xu — Xg) and {Xg,Xu) are independent. Therefore, from (|4.5p we 
get 

E{Xt\Tg,u) ==Xg + {Xu - Xg)^ 

u — s 

and thus X is a harness. The second formula in ()4.5p gives 
(4.6) Y^v{Xt\Tg,u)^{Xu-X^' {t-s){u-t) 



{U - S)2(c(u - S) + 1) 

Define now 



(4.7) Yt = ^c+^[{V + t)Xn^-t), te[0,oo). 

It is elementary to check that (^t)te[o,oo) is a quadratic harness and the param- 
eters are as follows 

2VV -2 

VT+W VvVT+cV 
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V 1 2 

^~ i + cv ' ~ v{i + cV) ' ~ (1 + cvy 



Var(7V,|^,,„) = ' " .^ ' {Nu N^). 



(Note that this agrees with the answers deduced from Proposition l4. II which imphes 
that Oy = 4ry, ry^ = Aay and 7y = 1 — 2^(Tyry.) 

On the other hand, it can be easily seen that the process X is a bridge on the 
gamma process (G't)te[o,oo) governed by the gamma distribution with the shape 
parameter 1/c and the scale equal 1. More precisely, in distribution process X 

is identical to the gamma bridge (Gt/Gv)te[o,v]\Gv = (Gt/Gv)tg[o,v], see [131 
Definition 2], see also [llj . 

Example 4.2 (Binomial process). For any real 7jo,6o, such that 770^0 = —l/N < 
for some iV g IN there exist a quadratic harness Y (namely, the Binomial process 
described here) which has parameters cry — Ty — 0, 9y = Oq, rjy — rjOi and 7y = 1. 

Indeed, consider conditional harnesses obtained as bridges from a Poisson pro- 
cess. Poisson process Nt with parameter A > is a harness with mean E,{Nt) — At 
variance Var(iVt) = Xt, and with conditional variance with respect to natural a- 
fields given by 

{u - sy 

Then 

Zt = Nt/^ - t 

is in gi?(0, 1;0,0;1). 

Instead of considering a conditional process of {Nt)t>o we therefore consider 
a conditional process of Z and apply Proposition 14.11 From (|4.2I) we see that 
7x = 1, crx = Tx =0. Indeed, to simplify the notation, consider i? = . Since 
Aov^ = (iV - V)/V for some e IN, from we see that Oy = -i^y = 1/Vn. 
So rjyOy = —1/N for some integer N E'iN. 

The conditional processes of the Poisson process are the Binomial processes. So 
the same conclusion can be obtained more directly without invoking explicitly any 
of the transformations. (Compare [51 Proposition 4.4].) Let b{n,p) denote the 
binomial distribution with sample size n and probability of success p. For fixed 

€ IN, define a Markov process X — {Xt)t£io,v] by the following (consistent) 
family of marginal and conditional distributions: 

Xt^b(^N,-^^ and Xt-Xs\Ys^b(^N-Ys,-^^^ , 0<s<t<V . 

Then the process X is called a binomial process with parameter N. It is elementary 
to see that the conditional distribution Xt\J-s,u ^ b (Yu — Kj, ■ Therefore X 
is a harness i.e. holds, and for any s,t,u E [0,V], s < t < u 

Var(X,| J-,.„) ^ l!ifM_i)(x„ - X,) . 

[u — s)^ 

An easy computation (or an application of Theorem II. ip shows that if 

{V + t)X_n_ ~tN 

Yt = , te[0, c»). 



then the process {Yt)t>Q is a quadratic harness and the parameters are 9y — y^V/N, 
rjy = —1/^/VN, Ty = ay ^ and 7y = 1. 
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On the other hand, it is immediate that X is a bridge obtained by conditioning 
a Poisson process (iVt)t>o, that is X = {Nt)teloy]\Nv = N. 

It is interesting to see which properties of a quadratic harness are preserved by 
conditioning. The only universal invariant that we found is related to parameter q 
of the Askey- Wilson law from the construction in fl] . 

Proposition 4.2. Suppose that Y e QH{rjY ^ Oy', <Jy , ty', 7y) is o, quadratic harness 
arising from a conditional bridge X m Z G QH(j]z, Qz'i<^z,tz',^z)- 

(i) 7z = ~1 */ o,nd only if "/y = ^1- 

(ii) //7z > —1, then 

(1 - 7y)^ - 4cryTy _ (1 - 7z)^ " "^CTZTz 

^ ' (l + 7y)2 " (l + 7z)2 

(iii) If Z is a q-Meixner process, i.e. rjz = (Jz = 0, "/ — q Cz [—1, 1] and the 
conditioning is with R = and Z]^ — 0, then 

(4.9) sign(4 - 4ty) - sign(02 _ 4^^) 

Proof. The first statement follows from (|1.23l) . Formula (14. 8p follows by a direct 
computation from (|1.21f[TT23)) . Formula (|4.9p follows from the expressions in Re- 
markO which gives 9^ - 4ty = (6*1 - 4tz)/(F + tz). □ 

Remark 4.1. For a Meixner process Z, the conclusion of Proposition iii) can be 
strengthened to double equality 

sign(7y^ - 4o-y) = sign(6'^ - 4ty) = sign(6'| - 4tz) . 



Indeed, from Proposition 14. II we see that — 4cry = (0y — 4:TY)crY /ty. 

In particular, quadratic harness Y that arises by conditioning from the negative 
binomial process Z has parameters riY^/W + Sy^/oy = 0, 7y = 1 — 2^cryTr, and 
d\ > Aty- (In fact, Y rises from a discrete "negative hypergeometric" processes.) 
Similarly, quadratic harness Y that arises from a conditional bridge of the "hyper- 
bolic Meixner" process Z with 0| < Atz has parameters Tyy-y/ry + BY\f^ = 0, 
7y = 1 — 2y/aYTY, and Oy < 4Ty. 

4.2. Proof of Theorem 11.31 From ()1.14p it is enough to construct a quadratic 
harness with parameters ax — tx & (0, 1), 9x — —itx and jx = 1 — 2crx. To do so, 
we determine the parameters of Z from (|4.2p and (14.31) . noting that for a Meixner 
processes A^y = {Zy — Zii)/{V — R) can take any positive real value that we fix 
first. 

5. Gluing construction 

This section is motivated by the construction of a classical bi-Poisson process 
from a pair of two conditionally independent Poisson processes |^ and by an- 
other recent construction of a quadratic harness from two conditionally indepen- 
dent copies of a negative binomial process [TB]. These constructions essentially 
consist of choosing an appropriate deterministic moment of time V so that Zy- 
conditional (and yy-conditionally independent) processes X+ := {Zt : t > V}\Zv 
and X_ := {Zt : t < V}\Zv arise as space time transforms of the above mentioned 
Levy processes. 



20 



WLODEK BRYC AND JACEK WESOLOWSKI 



We remark that in principle all quadratic harnesses arise from such a gluing 
construction. A fixed V > can be treated as the upper value V in Remark 
11.21 resulting in the process X_ = {Zt)t<v conditioned with respect to Zy, and 
on the other hand, we can treat this value as the left-hand side value R in Re- 
mark [L2] and consider the process X-)_ = {Zt)t>R conditioned with respect to Zj^. 
Then both processes are quadratic harnesses and by Markov property they are Zy- 
conditionally independent. So for example, a Poisson process arises from gluing 
a binomial process with another Poisson process, or a Wiener process arises from 
gluing a Brownian bridge with another Wiener process. 

The question of interest here is when the "components" X_ and X+ to be glued 
are in some sense "simpler" than the resulting process. Using Corollarv ll.2'1 and 
Remark 11.21 we can recognize when the components of such a gluing construction 
come from a well understood class of g-Meixner processes. (Since we are proving 
necessity only, we do not analyze X-|_.) 

Proposition 5.1. Let Z G QH{ri, 9; cr, r, 7) be defined on (0, od). Suppose that there 
is a deterministic moment of time V > such that Markov process X_ obtained as 
{Zt : t < V} conditioned with respect to Zy , can he transformed via (jl.27p into a 
q-Meixner process Y, then one of the following cases must happen: 

(i) 7 = 1,(T = T = and rj = 9 = 0. (Then X is the Wiener processes and the 
construction indeed works with any V > 0.) 

(ii) 7 = 1, fT = T = and rjO > 0. (Then V — 9 jr], X is the Poisson processes 
with parameter A which depends on Yy and the construction indeed works, 
see W .) 

(iii) fy-y/r — 9^fa, and 7 = 1 + 1\faT > 1. (Then V = \/t j o and X G 
Q-ff (0, ^x; 0, Tx; 1) with the sign of 9\ — Atx determined by the sign of 
9^ — At of process Y.j 

Proof. If Y comes from gluing a g-Meixner process then the conditional bridge Y_ 
corresponding to i? = 0, ^ = exists and can be transformed into a g-Meixner 
process X with parameters given in Remark ll.2f i). 

The only possibility for (|1.3ip to correspond to a q-Meixner process is when the 
parameters of Y satisfy 

(5.1) aV^ + {1--/)V + T = 0. 

Since cr, r > (see |3l Theorem 2.2]), the only solution with 7 < 1 is 7 = 1, 
a = T = 0. Then from (|1.29p we see that X is indeed a Mcixner process when we 
set V — 9/ri (which gives the bi- Poisson process) or when V" > but rj = 9 = 
(which gives the Wiener process). 

Other solutions of (|5.ip are possible only when (1 — 7)^ > 4(Tr. However, since 
7 < 1 -I- by [3l Theorem 2.2], this gives 7=1 + and V = ^/r/a. Then 

from (|1.29p . the coefficient at zy vanishes so X is indeed a Meixner process when 
ri^/r = 9^. 

For V — yjr/a^ 7 = 1 + formulas from Remark 11.21 (i) give 

e% _ + 46>Vg?zv + AcjTzl 
Atx ~ At + A9y/aTZy + Autz'^ ' 



so the sign of 9^ — At is preserved. 



□ 
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We remark that [16j provided gluing construction of two conditionally-independent 
copies of negative binomial processes. The resulting quadratic harness on (0, oo) 
reaches the "upper limit" 7=1 + l^/ar of the bound in [31 Theorem 2.2] and that 
for this process the product ctt can take arbitrarily values in (0, 00). We do not yet 
know whether processes in Proposition 15. II indeed exist for all signs of 9^ — 4t. 

Acknowledgement. We would like to thank Arthur Krener and Ofcr Zcitouni for 
information on reciprocal processes. Maja Jamiolkowska worked out the construc- 
tion [16J which inspired our Proposition 15.11 We also benefited from discussions 
with Ryszard Szwarc. This research was partially supported by NSF grant #DMS- 
0904720, and by Taft Research Seminar 2008/09. 

References 

[1] Wlodek Bryc and Jacek Wesolowski. Askey- Wilson polynomials, quadratic harnesses and 

martingales. Annals of Probability, 38:1221-1262, 2010. arxiv.org/abs/0812.0657. 
[2] Wlodzimierz Bryc, Wojciech Matysiak, Ryszard Szwarc, and Jacek Wesolowski. Projection 

formulas for orthogonal polynomials. arxiv.org/abs/math.CA/0606092, 2006. 
[3] Wlodzimierz Bryc, Wojciech Matysiak, and Jacek Wesolowski. Quadratic harnesses, q- 

commutations, and orthogonal martingale polynomials. Trans. Amer. Math. Soc, 359:5449- 

5483, 2007. arxiv.org/abs/math.PR/0504194. 
[4] Wlodzimierz Bryc, Wojciech Matysiak, and Jacek Wesolowski. The bi-Poisson process: a 

quadratic harness. Ann. Probab., 36:623-646, 2008. arxiv.org/abs/math.PR/0510208. 
[5] Wlodzimierz Bryc and Jacek Wesolowski. Conditional moments of q-Meixner processes. 

Probab. Theory Related Fields, 131:415-441, 2005. arxiv.org/abs/math.PR/0403016. 
[6] Wlodzimierz Bryc and Jacek Wesolowski. Classical bi-Poisson process: an invertible quadratic 

harness. Statist. Probab. Lett., 76:1664-1674, 2006. arxiv.org/abs/math.PR/0508383. 
[7] Wlodzimierz Bryc and Jacek Wesolowski. Bi-Poisson process. Infin. Dimens. Anal. Quantum 

Probab. Relat. Top., 10:277-291, 2007. arxiv.org/abs/math.PR/0404241. 
[8] Q.L. Burrell. burrelll988predictive. In L. EGGHE and R. ROUSSEAU, editors, Informet- 

rics 87/88: Select Proceedings of the First International Conference on Bibliometrics and 

Theoretical Aspects of Information Retrieval. Elsevier, 1988. 
[9] Q.L. Burrell. Modelling the Bradford phenomenon. Journal of Documentation, 44(1):1-18, 

1988. 

[10] Kjell Doksum. Tailfree and neutral random probabilities and their posterior distributions. 

The Annals of Probability, 2(2):183-201, 1974. 
[11] Michel Emery and Marc Yor. A parallel between Brownian bridges and gamma bridges. Publ. 

Res. Inst. Math. Sci., 40(3):669-688, 2004. 
[12] Thomas S. Ferguson. A bayesian analysis of some nonparametric problems. The Annals of 

Statistics, 1(2): 209-230, 1973. 
[13] Thomas S. Ferguson. Prior distributions on spaces of probability measures. The Annals of 

Statistics, 2(4):615-629, 1974. 
[14] J. M. Hammersley. Harnesses. In Proc. Fifth Berkeley Sympos. Mathematical Statistics and 

Probability (Berkeley, Calif., 1965/66), Vol. Ill: Physical Sciences, pages 89—117. Univ. 

California Press, Berkeley, Calif., 1967. 
[15] Jean Jacod and Philip Protter. Time reversal on Levy processes. Ann. Probab., 16(2):620-641, 

1988. 

[16] Maja Jamiolkowska. Bi-Pascal process - definition and properties. Master's thesis, Warsaw 

University of Technology, (in Polish) 2009. 
[17] B. Jamison. Reciprocal processes. Probability Theory and Related Fields, 30(l):65-86, 1974. 
[18] Roger Mansuy and Marc Yor. Harnesses, Levy bridges and Monsieur Jourdain. Stochastic 

Process. Appl, 115(2):329-338, 2005. 
[19] Jacek Wesolowski. Stochastic processes with linear conditional expectation and quadratic 

conditional variance. Probab. Math. Statist., 14:33-44, 1993. 
[20] David Williams. Some basic theorems on harnesses. In Stochastic analysis (a tribute to the 

memory of Rollo Davidson), pages 349-363. Wiley, London, 1973. 



22 



WLODEK BRYC AND JACEK WESOLOWSKI 



[21] Evdokia Xckalaki and Mimoza Zografi. The generalized Waring process and its application. 

Comm. Statist. Theory Methods, 37(11-12);1835-1854, 2008. 
[22] M. Zografi and E. Xckalaki. The generalized Waring process. In Proceedings of the 5th 

Hellenic-European Conference on Computer Mathematics and its Applications, Athens, 

Greece, pages 886-893, 2001. 

Appendix A. Auxiliary details 

Department of Mathematics, University of Cincinnati, PO Box 210025, Cincinnati, 
OH 45221-0025, USA 

E-mail address: Wlodzimierz.BrycaUC.edu 

Faculty of M;\thematic;s and Information Science, Warsaw University of Technol- 
ogy, PL. POLITECHNIKI 1, 00-661 WARSZAWA, POLAND 
E-mail address: wesoloSalpha. mini. pw.edu. pi 



